Jet measurements at D0 using a algorithm 
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D0 has implemented and calibrated a fc^ jet algorithm for the first time in a pp collider. We present two 
results based on 1992-1996 data which were recently published; the subjet multiplicity in quark and gluon jets 
and the central inclusive jet cross section. The measured ratio between subjet multiplicities in gluon and quark 
jets is consistent with theoretical predictions and previous experimental values. NLO pQCD predictions of the 
fcx inclusive jet cross section agree with the D0 measurement, although marginally in the low pT range. We also 
present a preliminary measurement of thrust cross sections, which indicates the need to include higher than 
terms and resumation in the theoretical calculations. 



1. INTRODUCTION 

Until recently, only cone algorithms were used 
to reconstruct jets in pp [Q colliders. The cone al- 
gorithm used to reconstruct the 1992-1996 Teva- 
tron data presents several shortcomings: 

• An arbitrary procedure must be imple- 
mented to split and merge overlapping 
calorimeter cones. 

• An ad-hoc parameter, i?sGp @|, is required 
to accommodate the differences between jet 
definitions at the parton and detector levels. 

• Improved theoretical predictions calcu- 
lated at the next-to-next-to-leading-order 
(NNLO) in pQCD are not infrared safe, but 
exhibit sensitivity to soft radiation. 

A second class of jet algorithms, free of 
these problems, was developed during the past 
decade These recombination algorithms 

successively merge pairs of nearby objects (par- 
tons, particles, or calorimeter towers) in order of 
increasing relative transverse momentum. A sin- 
gle parameter, Z?, which approximately charac- 
terizes the size of the resulting jets, determines 
when this merging stops. No further splitting 
or merging is involved because each object is 
uniquely assigned to a jet. There is no need to in- 
troduce any ad-hoc parameter, because the same 



algorithm is applied at the theoretical and exper- 
imental level. Furthermore, by design, clustering 
algorithms are infrared and coUinear safe to all 
orders of calculation. 

The D0 Collaboration has implemented a fcx 
algorithm to reconstruct jets from data taken dur- 
ing the 1992-1996 collider run. This paper is a 
review of the associated measurements recently 
performed by the D0 experiment |8[|^ . 



2. THE RUN I D0 DETECTOR 

D0 is a multipurpose detector designed to 
study pp collisions at the Fermilab Tevatron Col- 
lider. A full description of the Run I D0 detector 
can be found in Ref. The primary detector 

components for jet measurements at D0 are the 
calorimeters, which use liquid-argon as the ac- 
tive medium and uranium as the absorber. The 
D0 calorimeters provide full solid angle coverage 
and particle containment (except for neutrinos or 
high pt muons), as well as linearity of response 
with energy and compensation (b/tt response ra- 
tio is less than 1.05). Transverse segmentation is 
Ar]xA(f> = 0.1x0.1, with?7 = -In [tan(6'/2)]. The 
single particle energy resolutions for electrons (e) 
and pions (tt) , measured from test beam data, are 
approximately 15% and 50% respectively. 
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3. fc_L JET ALGORITHM 

The D0 k± jet algorithm |^ starts with a Ust 
of energy pre-clusters, formed from calorimeter 
cells, final state partons or particles. The an- 
gular separation between pre-clusters is ATZ = 
^/A^f + A^ > 0.2. 



(a) (b) 




* 



Figure 1. Example of the k± merging mechanism 
on a hard pp event, (a) The particles represented 
by arrows comprise a list of objects, (b-f) Solid 
arrows represent the final jets reconstructed by 
the k± algorithm, and open arrows represent ob- 
jects not yet assigned to jets. The five diagrams 
show successive iterations of the algorithm. In 
each diagram, an object is labeled as a jet if it 
is well separated from all other objects. Two ob- 
jects are merged if they have small relative k±. 
The asterisk labels the relevant objects at each 
step. 



The jet reconstruction is performed in three 
steps: 

1. For each object i in the list, the algorithm 
defines = p^j, where p^ is the momentum 
transverse to the beam. For each pair of 

objects, it also defines = min(p^;,p^j) ^^'^ , 



where D is a resolution parameter. 

2. If the minimum of all possible dii and dij is 
a dij, the algorithm replaces objects i and j by 
their 4- vector sum before going back to step 1. If 
the minimum is a da, then i is defined as a jet 
and removed from the list of objects. 

3. While there is an object left in the list, the 
algorithm returns to step 1. 

The final product of this process is a list of 
jets, separated by ATZ > D from each other. 
Subjets may be defined by re-running the k± 
algorithm from a list of pre-clusters in a given 
jet. Pairs of objects with the smallest are 
merged successively until all remaining dy satisfy 
dij > ycutPxjcf The resolved objects are called 
subjets, and the number of subjets within the jet 
is the subjet multiplicity M. For j/cut==l, the en- 
tire jet consists of a single subjet (M=l). As j/cut 
decreases, the subjet multiplicity increases, until 
every pre-cluster becomes resolved as a separate 
subjet in the limit ycut ^ 0. 

4. CALIBRATION OF JET MOMEN- 
TUM 

The uncertainty in the jet energy or momen- 
tum is the dominant error in almost every jet 
measurement at a hadron collider. The jet mo- 
mentum calibration is described in Ref. 
The calibration at D0 accounts for detector ef- 
fects like response, noise, and signal pile-up from 
previous crossings. It also removes the under- 
lying event formed by the remnant soft partons 
(u.e.), and the contribution of multiple pp interac- 
tions. These corrections enter a relation between 
the momentum of a jet measured in the calorime- 
ter p^^'^^ and the "true" jet momentum Pjg"°: 

t.ue J^]gr-Po(V"^Apjr") 

i?jet(77j<=\p^<=t) ^ ' 

where the offset term po corrects for u.e., noise, 
pile-up, and multiple interactions, while i?jot cor- 
rects for the response of the calorimeter to jets. 
The true jet momentum is defined as the parti- 
cle level jet momentum. A particle level jet is 
reconstructed from the final state particles, af- 
ter hadronization but before interaction with the 
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calorimeter material. The calibration procedure 
follows closely that of the calibration of the fixed- 
cone jet algorithm ||lT[| . 

The fractional momentum resolution for k± 
jets {D=l) is determined from the mea- 
sured pt imbalance in dijet events. At 
100 (400) GeV, the fractional resolution is 
0.061±0.006(0.039±0.003). Within statistical 
and systematic uncertainties, there is not a sig- 
nificant difference between energy resolutions as- 
sociated with k± {D=l) and cone TZ=0.7 jets. 

5. PHYSICS RESULTS 

D0 has performed a number of measurements 
using the k± jet algorithm. These include a study 
of the structure of quark and gluon jets and 
measurements of the central {\ri\ < 0.5) inclusive 
jet cross section and thrust distributions. 

5.1. Subjet multiplicities 

In LO QCD, the fraction of final state jets 
which are gluons decreases with x Pt/Vs, 
the momentum fraction of initial state partons 
within the proton. For fixed px, the gluon jet 
fraction decreases when ^/s is decreased from 
1800 GeV to 630 GeV. We select gluon and quark 
enriched samples with identical cuts in events at 
-^i = 1800 and 630 GeV to reduce experimental 
biases and systematic effects. Of the two highest 
Pt jets in the event, we select k±{D = 0.5) jets 
with 55 < PT < 100 GeV and |ry| < 0.5. 

There is a simple method to extract a measure- 
ment of quark and gluon jets on a statistical basis. 
If M is the subjet multiplicity in a mixed sample 
of quark and gluon jets, it may be written as a 
linear combination of subjet multiplicity in gluon 
and quark jets: 



M - /7\fg + (1 - /)Mq 



(2) 



The coefficients are the fractions of gluon and 
quark jets in the sample, / and (1 — /), respec- 
tively. Consider Eq. (||) for two similar samples 
of jets at ^/s — 1800 and 630 GeV, assuming Afg 
and Mq are independent of y^. The solutions are 



Mg = 



(l - /*^30) il/lSOO - (1 - /18"0) 1/63" 



^1800 _ f 



630 



(4) 



where M^^°" and M^^o g^j.g ^j^g experimen- 
tal measurements in the mixed jet samples at 
= 1800 and 630 GeV, and f^^°° and are 
the gluon jet fractions in the two samples. The 
method relies on knowledge of the two gluon jet 
fractions, which are extracted from the herwig 
5.9||lJ Monte Carlo event generator and used in 
Eqs^lfV 

Figurcl2 shows the average subjet multiplicity 
(z/cut = 10^3) for quark and gluon jets. Mg is sig- 
nificantly larger for gluon jets than for quark jets. 
The gluon jet fractions are the dominant source 
of systematic error. We also compute the ratio 



p _ (M )-l _ 



= 1.84±0.15(stat)+°;^g(sys). Fig- 



ure m shows a comparison between the ratio mea- 
sured by D0 , the herwig 5.9 result of r=1.91, 
the ALEPHII value of r=1.7±0.1 (e+e" anni- 
hilations at \fMz — Mz) , and the associated 
Monte Carlo and resummation prediction [p^ . 
Good agreement is observed. All of the exper- 
imental and theoretical values for r are smaller 
than the naive QCD prediction value of 2.25 for 
the ratio of color charges. This is because of 
higher-order radiation in QCD, which tends to 
reduce the ratio from the naive value. 
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Figure 2. Corrected subjet multiplicity for gluon 
and quark jets, extracted from D0 data. 
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Figure 3. The subjet multiplicity in (a) gluon 
and (b) quark jets, for D0 data, for the herwig 
Monte Carlo, and resummed predictions. 

5.2. Central inclusive jet cross section 

The inclusive jet cross section for |7/| < 0.5, 
d'^a/{dpTdri), was measured as N/{AriApT€L), 
where A77 and Apx are the rj and pt bin sizes, 
N is the number of jets reconstructed with the 
k± (D—l) algorithm in that bin, e is the over- 
all efficiency for jet and event selection, and L 
represents the integrated luminosity of the data 
sample 

The fully corrected cross section for |r/| < 0.5 
is shown in Fig. ^, along with the statistical un- 
certainties. The systematic uncertainties include 
contributions from the jet and event selection, un- 
smearing, luminosity, and the uncertainty in the 
momentum scale, which dominates at all trans- 
verse momenta. The fractional uncertainties for 
the different components are plotted in Fig. || as 
a function of the jet transverse momentum. 

The results are compared to the NLO pQCD 
prediction from jetrad [0, with the renormal- 
ization and factorization scales set to p^^^/2, 
where p^'^^ refers to the px of the leading jet in 
an event. The comparisons are made using pa- 
rameterizations of the parton distribution func- 
tions (PDFs) of the CTEQ H and MRST 
families. Figure || shows the ratios of (data- 
theory)/theory. The predictions lie below the 
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Figure 4. The central {\ri\ < 0.5) inclusive jet 
cross section obtained with the k± algorithm at 
y/s = 1.8 TeV. Only statistical errors are in- 
cluded. The solid line shows a prediction from 
NLO pQCD. 



data by about 50% at the lowest pT and by (10- 
20)% for PT > 200 GeV. To quantify the compari- 
son in Fig. |6| the fractional systematic uncertain- 
ties are multiplied by the predicted cross section, 
and a comparison, using the full correlation 
matrix, is carried out Though the agreement 
is reasonable (x^/dof ranges from 1.56 to 1.12, 
the probabilities from 4 to 31%), the differences 
in normalization and shape, especially at low px, 
are quite large. The points at low px have the 
highest impact on the x^- If the first four data 
points are not used in the comparison, the 
probability increases from 29% to 77% when us- 
ing the CTEQ4HJ PDF. 
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the cross section. 
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Figure 6. Difference between data and JE- 
TRAD pQCD, normalized to the predictions. 
The shaded bands represent the total system- 
atic uncertainty. In the bottom plot a herwig 
hadronization contribution has been added to the 
prediction (open circles). 



The NLO predictions of the inclusive cross 
section for kj_ {D=l) and cone jets (7?.=0.7, 
i?sep=l-3) in the same |r/| < 0.5 interval are 
within a few percent of each other in the pt 
range relevant in this analysis The mea- 

sured k± cross section, however, is 37% (16%) 
higher than the published cone algorithm ||2^ at 
60 (200) GeV. This difference in the cross sections 
is consistent with the measured difference in px 
for cone jets matched in r]-(j) space to k± jets, as 
shown in Fig. ^. For the same energy clusters 
the Pt of k± jets (D—l) is higher than the Et 
of associated cone jets {TZ=0.7). The difference 
increases with jet pt, from about 5 GeV (or 5%) 
at PT « 100 GeV to about 7 GeV (or 3%) at 
Pt « 250 GeV ||]. Fig. ^ proves, however, that 
the energy difference does not depend on the in- 
stantaneous luminosity associated with the sam- 
ple. After offset subtraction, it is clear that k± 
jets are not contaminated by energy coming from 
pile-up, uranium noise, or multiple interactions. 

The effect of final-state hadronization on re- 
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Figure 7. pt of k± jets {D=l) minus Et of the 
associated cone jets (7?.=0.7) for samples taken at 
different instantaneous luminosities. 

constructed energy, which might account for the 
discrepancy between the observed cross section 
using k± and the NLO predictions at low pT, 
and also for the difference between the k± and 
cone results, was studied using herwig (version 
5.9) simulations. Figure js] shows the ratio of pt 
spectra for particle-level to parton-level jets, for 
both the k± and cone algorithms. Particle cone 
jets, reconstructed from final state particles (af- 
ter hadronization), have less pt than the parton 
jets (before hadronization), because of energy loss 
outside the cone. 
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Figure 8. Ratio of particle-level over parton-level 
HERWIG PT spectra for jets, as a function of the 
parton jet transverse momentum. 

In contrast, k±^ particle jets are more energetic 
than their progenitors at the parton level, due to 
the merging of nearby partons into a single par- 
ticle jet. Including the hadronization effect de- 
rived from HERWIG in the NLO jetrad predic- 
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tion improves the probability from 29% to 44% 
(31% to 46%) when using the CTEQ4HJ (MRST) 
PDF. 

We have also investigated the sensitivity of the 
measurement to the modeling of the background 
from spectator partons through the use of mini- 
mum bias events, and found that it has a small 
effect on the cross section: at low px, where the 
sensitivity is the largest, an increase of as much as 
50% in the underlying event correction decreases 
the cross section by less than 6%. 

5.3. Thrust cross sections 

Event shape variables have been extensively 
used in e~ and ep collider experiments to study 
the spatial distribution of hadronic final states, to 
test the predictions of perturbative QCD, and to 
extract a precise value of the coupling constant 
Qfs- Over the last few years, they have attracted 
considerable interest, as they have proved to be 
a fruitful testing ground for recent QCD devel- 
opments like resummation calculations and non- 
perturbative corrections. 

There are several observables which character- 
ize the shape of an event. To be calculable by per- 
turbation theory, these quantities must be infra- 
red safe, i.e. insensitive to the emission of soft 
or collinear gluons. A widely used variable that 
meets this requirement is the thrust, defined as 



T 



maxft- 



(5) 



where the sum is over all partons, particles or 
calorimeter towers in the event. The unit vector 
n that maximizes the ratio of the sums is called 
the thrust axis. The values of thrust range from 
T=0.5 for a perfectly spherical event, to T=l for 
a pencil-like event, when all emitted particles are 
collinear. In this latter case, the thrust axis lies 
along the direction of the particles. 

In most of the kinematic range, e^e~ and ep 
collider experiments ||2l|] report good agreement 
of event shape distributions with 0[a^) pQCD 
corrections to the lowest order QED diagram that 
governs the interaction. Fixed order QCD calcu- 
lations, however, fail when two widely different 
energy scales are involved in the event, leading to 
the appearance of large logarithmic terms at all 



orders in the perturbative expansion [g2j. This 
happens in the limit of the two jet back-to-back 
configuration, when T ^ 1. This case is handled 
by the resummation technique, which identifies 
the large logarithms in each order of perturba- 
tion theory and sums their contributions to all 
orders. DELPHI reports excellent agreement of 
thrust distributions in Z ^ hadrons once resum- 
mation and hadronization corrections are added 
to the 0(a2) QCD prediction ||. 

In a hadron collider, it is convenient to intro- 
duce "transverse thrust", T"^, a Lorenz invariant 
quantity under z-boosts, which is obtained from 
Eq. H in terms of transverse momenta: 

Ei \PT, ■ n\ 



T = maxfi- 



Ei Ipt.I 



(6) 



Transverse thrust ranges from T'^=l to 
r"^=2/7r ((I cos 6'!)) for a back-to-back and an 
isotropic distribution of particles in the transverse 
plane, respectively. To minimize systematics as- 
sociated with the busy environment in a pp col- 
lider, we use only the two leading jets in the event, 
reconstructed with a k± D—1 algorithm, rather 
than using all calorimeter towers. Other particles 
in the event are inferred from the angular distri- 
bution of the two leading jets. 

The measurement of the dijet transverse thrust, 
T2 presents a good opportunity to test resum- 
mation models, as well as the recently developed 
NLO pQCD three-jet generators ^]^. is 
binned in terms of i?T3, defined as the scalar sum 
of the transverse momenta of the three leading 
jets of the event pxi +PT2 +PT3- Hts is an esti- 
mator of the energy scale of the event. The lowest 
order at which pQCD does not give the trivial re- 
sult T2—I is 0{a^), corresponding to up to three 
parton jets in the final state. A 0{a^) calcula- 
tion, like JETRAD, does not cover the whole phys- 
ical range of T2 . A 0{a^) prediction will also fail 
at T2 1, where soft radiation in dijet events 
contribute large logarithms which need to be re- 
summed. 

Figures ^ and |lO| show thrust cross sections as 
a function of 1 — T]- for the lowest and highest 
Ht3 bins: 160-260 GeV and 430-700 GeV. 

The error bars are statistical, and the system- 
atic uncertainties are in the 15-30% range. The 
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Figure 9. Thrust cross section as a function of 
1 - T]^ for the lowest Ht3 bin: 160-260 GeV. 
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Figure 10. Thrust cross section as a function of 
1 - T]^ for the highest Hts bin: 430-700 GeV. 



higher the value of 1 — T2 , the more dominant is 
the contribution of high order terms in as- This 
explains the disagreement between the data and 
the prediction jetrad in the high 1 — T2 
range. For example, for between \/2/2 and 
the Qfg terms contribute to LO. When we 
expand the range of thrust around unity using 
a logarithmic scale, we verify that jetrad also 
fails, suggesting the need to resum higher order 
terms to improve agreement. Table ^ displays the 
agreement probability between the D0 data and 
JETRAD, which decreases sharply as we incorpo- 
rate data points in the high and low I — T2 range. 
The probabilities are calculated using the full co- 
variance error matrix in a test. 



Table 1 

Agreement probability between the D0 data and 
JETRAD, using the full covariance error matrix in 



a X test. 
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6. Conclusions 

D0 has successfully implemented and cali- 
brated a k± jet algorithm in a pp collider. Quark 
and gluon jets have a different structure consis- 
tent with the HERWIG prediction and previous ex- 
perimental results from — e~ colliders. The 
thrust cross section measurements indicate the 
need to include higher than terms and resuma- 
tion in the theoretical predictions. They also offer 
an excellent opportunity to test the recently de- 
veloped NLO three jet generators |24,2^. The 
marginal agreement between the measurement of 
the particle level inclusive k± jet cross section 
with the ttg theory is opening a debate on mat- 
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ters such us hadronization, underlying event, and 
algorithm definition. 
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